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Abstract 

A free-piston Stirling engine performance 
code is being upgraded and validated at the NASA 
Lewis Research Center under an interagency agree- 
ment between the Department of Energy's Oak Ridge 
National Laboratory and NASA Lewis. Many modi- 
fications have been made to the free-piston code 
in an attempt to decrease the calibration effort. 

A procedure was developed that made the code cali- 
bration process more systematic. Engine-specific 
calibration parameters are often used to bring 
predictions and experimental data into better 
agreement. The code was calibrated to a matrix of 
six experimental data points. Predictions of the 
calibrated free-piston code are compared with 
RE-1000 free-piston Stirling engine sensitivity 
test data taken at NASA Lewis. Reasonable agree- 
ment was obtained between the code predictions 
and the experimental data over a wide range of 
engine operating conditions. 

Executive Summary 

A free-piston Stirling engine performance 
code is being upgraded and validated at NASA Lewis 
under an interagency agreement between the Depart- 
ment of Energy's Oak Ridge National Laboratory and 
NASA Lewis. Many modifications have been made to 
the free-piston code in an attempt to decrease the 
calibration effort required. The following changes 
have been made to the code: (1) heat-transfer and 

f riction-f actor correlations for the cooler and 
the regenerator were modified; (2) gas leak path 
model for the power piston was modified; (3) dis- 
placer appendix gap pumping loss equation was 
added; (4) option was added to separate the con- 
necting ducts into their own control volumes; and 
(5) option was added to run the code with the dis- 
placer gas spring represented by either a control 
volume or a spring constant and damping factor. 

After all the changes were made to the code, 
a process was developed that logically arrived at 
the best set of engine-specific calibration param- 
eters. Engine-specific calibration parameters 
based on engine test data were used to bring pre- 
dictions and experimental data into better agree- 
ment. These parameters were defined as the set of 
multiplication factors and coefficients required 
to adjust predicted pressure drops, heat transfer, 
and gas flow rates so that the code predictions 
could better agree with a specific engine's test 
data. The code was calibrated to a matrix of six 
experimental data points. Predictions of the cali- 
brated code were compared with RE-1000 free-piston 
Stirling engine experimental data. Measured and 
predicted data that were compared include indi- 
cated power, indicated efficiency, compression space 
pressure amplitude, displacer stroke, piston- 
displacer phase angle, piston-compression space 
pressure phase angle, frequency, and expansion and 
compression space gas temperatures. Most compar- 
isons agreed within the error bands of the test 
data. The predicted indicated power and efficiency 


agreed within ±3 to ±11 percent of the actual 
test data. 

More information is needed in the areas of 
oscillating flow heat transfer and pressure drop, 
leakage and centering port flow, and gas spring 
hysteresis loss. The ultimate goal is to have a 
code that does not require any engine-specific 
calibration parameters. Although the predictions 
of the code agreed fairly well with the RE-1000 
test data, it is an engine-specific code. Before 
the code could be used to simulate a different 
free-piston engine with a large degree of confid- 
ence, the code would need to be calibrated for 
that particular engine. 

Nomenclature 

Ad cross-sectional area of displacer, in. 2 
(cm 2 ) 

A dr cross-sectional area of displacer rod, 
in. 2 (cm 2 ) 

A e ff effective area for heat transfer on water 
side of cooler, in. 2 (cm 2 ) 

Ap cross-sectional area of power piston, 

in. 2 (cm 2 ) 

C conversion factor 

C^s damping factor of displacer gas spring, 

lbf-s/in. (N-s/cm) 

D displacer diameter, in. (cm) 

Dfo hydraulic diameter, in. (cm) 

F Fanning friction factor 

Fg body force on displacer, lbf (N) 

F,j net force acting on displacer, lbf (N) 

F s displacer gas spring force, lbf (N) 

H heat-transfer coefficient, Btu/ft 2 -s-°R 

(W/m 2 -°C) or in.-l bf/ in. 2 -s-°R (W/cm 2 -°C) 

h radial clearance between displacer and 

displacer cylinder, in. (cm) 

^ds spring constant of displacer gas spring, 
lbf/in. (N/cm) 

^dsal displacer gas spring adiabatic spring 

constant limit, lbf/in. (N/cm) 

k aq thermal conductivity of gas in appendix 

gap, Btu/in.-s-°R (W/cm- C) 

k a ] thermal conductivity of aluminum, 

Btu/ in.-s-°R (W/cm- C) 
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k r thermal conductivity of gas in regenerator, 
in.-lb-f/in.-s-^R (W/cm-°C) 

k ss thermal conductivity of stainless steel 
displacer, Btu/in.-s-°R (W/cm- C) 

k w thermal conductivity of water, Btu/ft-s-°R 

W <W/m-°C) 

L length of cooler tubes, in. (cm) 

M d mass of displacer, slugs (kg) 

aP range of expansion space pressure, psi 

(N/cm 2 ) 

Pg mean bounce space pressure, psi (N/cm 2 ) 

Pq mean compression space pressure, psi 

(N/cm 2 ) 

P ra compression space pressure amplitude, psi 
(N/cm 2 ) 

Pq mean displacer gas spring pressure, psi 

(N/cm 2 ) 

P d instantaneous displacer gas spring 

pressure, psi (N/cm 2 ) 

P^ mean expansion space pressure, psi 

(N/cm 2 ) 

Pr Prandtl number 

p regenerator porosity 

Q out heat rejected from cooler per engine 

cycle, in. -lbf /cycle (cm-N/cycle) 

Qpump appendix gap gas enthalpy transfer (or 
pumping loss), in.-lbf/s (cm-N/s) 

R a l thermal resistance for conductive heat 

transfer in aluminum, s-°R/Btu (°C/W) 

Re Reynolds number 

RH20 thermal resistance for convective heat 

transfer between water and cooler wall, 

s-°R/Btu (°C/W) 

r^ effective radius of gas side of cooler, 

in. (cm) 

r 0 effective radius of water side of cooler, 

in. (cm) 

aT temperature differential, °R (°C) 

TH20 cooling water inlet temperature, °R (°C) 

maximum temperature of gas entering the 
appendix gap, °R (°C) 

T] minimum temperature of gas entering the 
appendix gap, °R (°C) 

T w temperature of the cooler wall, °R (°C) 

T wn i a temperature of the cooler wall during the 

previous cycle, °R (°C) 


Vdo mean displacer gas spring volume, in. 3 
(cm 3 ) 

V D amplitude of the piston velocity, in./s 
(cm/s) 

X d amplitude of displacer stroke, in. (cm) 

x d instantaneous position of displacer, in. 
(cm) 

x d instantaneous velocity of displacer, in./s 

(cm/s) 

x d instantaneous acceleration of displacer, 

ft/s 2 (m/s 2 ) 

Xp amplitude of piston stroke, in. (cm) 

Xp instantaneous position of piston, in. (cm) 

a thermal diffusivity of displacer, in. 2 /s 

(cm 2 /s ) 

y specific heat ratio, C p /C v 

9 phase angle between the power piston and 

the compression space pressure wave, deg 

<$> phase angle between the displacer and the 

expansion space pressure wave, deg 

phase angle between the compression space 
pressure wave and the piston velocity, deg 

w engine frequency, Hz 

Introduction 

A kinematic Stirling engine performance code 
written at NASA Lewis 3 was modified in the areas 
of engine thermodynamics, piston/displacer dynam- 
ics, and engine generic geometry to predict the 
performance of a free-piston Stirling engine. 2 
Mechanical Technology, Inc., (MTI) modified the 
kinematic code in 1980 under contract with NASA 
Lewis. Development of the free-piston code con- 
tinued at NASA under an interagency agreement 
between NASA Lewis and the Department of Energy's 
Oak Ridge National Laboratory (ORNL). Under this 
agreement, the RE-1000 free-piston Stirling engine 
has been modeled and tested. 

The RE-1000 was designed and built by Sun- 
power Inc., of Athens, Ohio. It was designed as a 
research tool and was obtained in 1979 for testing 
as part of the NASA Stirling engine technology 
program at Lewis. The RE-1000 is an electrically 
heated single-cylinder engine with a dashpot load. 
A cutaway view of the engine is shown in Fig. 1. 
The engine was tested over a wide range of heater- 
tube outside-wall temperatures, mean working-space 
pressures, cool ing-water inlet temperatures, pis- 
ton strokes, and working fluids. In addition, 
tests were conducted with different engine config- 
urations. These different configurations included 
two different regenerators, two displacers, and 
two pistons. When the engine was purchased, it 
was built with a regenerator and a displacer that 
are referred to as regenerator 1 and displacer 1. 
These parts were optimized within a set of design 
constraints to achieve high engine efficiency. At 
a later time in the testing of the RE-1000, NASA 
Lewis purchased a regenerator and displacer from 
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1500 < Re < 20 000 


Sunpower that were optimized for high-power output 
with all other parameters of the engine held con- 
stant. These high-power parts are referred to as 
regenerator 2 and displacer 2. For this report 
the high efficiency engine configuration was chosen 
and the code was calibrated. Comparisons were 
made between predicted and measured data over a 
range of operating conditions. 3 Regenerator 1 
had a porosity of 75.9 percent while displacer 1 
was designed to operate with a phase angle of 
about 45 with respect to the power piston. The 
actual phase angle at engine design conditions 
was higher than the design phase angle, probably 
due to a slightly higher displacer damping than 
assumed during the design process. The standard 
piston mass gave the engine a resonant frequency 
of about 30 Hz. 

This report describes the modifications made 
to the code since the validation effort was car- 
ried out as described by Tew. In addition, 
the approach used to calibrate the code is pre- 
sented. Finally, the code predictions are com- 
pared with experimental data for the RE-1000 
free-piston Stirling engine. 

Code Modifications 


The free-piston Stirling engine code in its 
previous form could predict trends in engine per- 
formance. Predicted power and efficiency were in 
reasonably good agreement with test data within 
the ±10 to ±15 percent range. In an effort to 
improve the code predictions, many modifications 
were made. These modifications are described in 
detail in the following paragraphs. 

The first modification involved the cooler 
model. For the RE-1000, the cooler is a cylindri- 
cal aluminum heat exchanger which consists of 135 
rectangular gas passages and 80 circular water 
passages. The cooler is both a parallel and 
counterflow type of heat exchanger depending on 
the instantaneous direction of the oscillating 
working space gas. All the passages run in an 
axial direction. The gas passages are located 
near the inner radius of the cylinder while the 
water passages are located near the outer radius 
of the cylinder as shown in Fig. 2. Previously, 
the gas side of the cooler was modeled using cir- 
cular tube heat-transfer and friction-factor 
correlations. 

Since the RE-1000' s cooler passages have 
rectangular cross sections, rectangular tube heat- 
transfer and friction-factor correlations were 
substituted for the circular tube correlations. 
Both the heat-transfer and the friction-factor 
correlations assume fully developed, steady flow 
and were obtained from Kays and London. 5 The 
rectangular passage heat-transfer correlation is 
in tabular form for a range of Reynolds numbers 
from 0 to 30 000. The method of linear inter- 
polation is used for Reynolds numbers that fall 
between the table values. The friction factor 
correlations for the gas side of the cooler are 

F = Re < 1500 

F = Re > 20 000 

Re U *^ 


c 0.079 

The friction-factor equation for laminar flow 
(i.e.. Re < 1500) is for the RE-1000's rectangular 
cooler passages which have an aspect ratio of 0.135. 
This equation must be revised if any changes are 
made to the dimensions of the cooler gas passages. 
The term aspect ratio is defined as 

aspect ratio = b/a b \~ [ 

a 

The friction-factor equations for turbulent flow 
are valid for any rectangular passages. 

A heat-transfer correlation was added for the 
water side of the cooler. This correlation elim- 
inated the need to guess the initial temperature 
of the cooler wall as was done previously. The 
cooling-water inlet temperature and flow rate, 
which can be accurately measured, are now used as 
inputs. The heat-transfer correlation for the 
water side of the cooler is 

C*k 0.0668(D,/L)RePr 

H = -R-* 3.66+ - n 

D h [ 1 + 0.04((D h /L)RePr) U - bbb 

where C = 12 in. /ft or 100 cm/m. This correlation 
assumes fully developed laminar flow in tubes with 
constant wall temperatures. 6 The heat-transfer 
coefficient calculated from this equation is then 
used in the calculation of the thermal resistance 
of the water side of the cooler as shown below. 
Also, the thermal resistance of the aluminum sep- 
arating the water and the gas is calculated assum- 
ing an effective thickness of material as shown 
below. The effective area and the effective thick- 
ness used in the thermal resistance equations are 
calculated assuming solid cylinderical water and 
gas passages with the same flow areas of the actual 
water and gas passages: 

R H2o = c Ha^: 

where C = 144 in.^/ft^ or 10 000 and 

1 n ( r a / r i ) 

K al “ 2 " Lk al 

After the thermal resistances are calculated, the 
temperature of the cooler wall on the gas side is 
calculated using the equations 

T w = T H20 + Q out (1>C ^ R al + R H20^ 


where C = 1.0725xl0“ 4 Btu/in.-lbf or 0.01 m/cm and 


T = 0.5T + 0.5T 

w w wold 


A correction was made to the code concerning 
the leak path for the power piston. The previous 
leakage model assumed that the major leak path for 
the piston was from the centering port on the 
piston to the buffer space through the annular 
clearance between the piston and the cylinder as 
shown in Fig. 3. The new leakage model simulates 
the leak path between the centering port on the 
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piston and the centering port on the piston cylin- 
der. The new leak path, as shown in Fig. 3, is 
much shorter than the old path. 


A displacer appendix gap gas enthalpy trans- 
fer (or pumping loss) equation has been added to 
the free-piston code. The term appendix gap 
refers to the small clearance space between the 
displacer and the displacer cylinder. The net 
appendix gap loss is now calculated by adding the 
net gas enthalpy loss to the shuttle loss. The 
shuttle loss is a conduction loss enhanced by the 
displacer's oscillatory motion and the radiation 
heat transfer between the displacer and cylinder 
walls; the net gas enthalpy loss is the heat 
transfer down the gap by virtue of the working 
gas motion, pressure, and temperature. This is 
illustrated in Fig. 4. The equation used for cal- 
culating the net gas enthalpy transfer was taken 
from Ref. 7 and is as follows: 


Q 


pump 


= irDAPhX^w 



In 




n”5 

h(u)/2a) U *\ 


I sin <t> 


The maximum temperature of the gas entering the 
appendix gap T^ was assumed to be at the 
cylinder hot end metal temperature while the min- 
imum temperature T] was assumed to be at the 
cylinder cold end temperature. This equation 
neglects the hysteresis loss that occurs in the 
appendix gap. In compensation, the surface area 
of the displacer cylinder in the appendix gap was 
included in the expansion space surface area for 
heat-transfer calculation. 


The working space model consists of the swept 
volumes in the expansion and compression spaces, 
three heat exchangers (heater, regenerator, and 
cooler), and four connecting ducts. The heater, 
regenerator, and cooler are connected in series 
via connecting ducts between the two swept volumes. 
The working space is divided into a number of con- 
trol volumes for analysis of fluid flow and heat 
transfer. Previously, the connecting duct volumes 
were lumped with adjacent heat exchangers or swept 
volumes. An option has been added to separate 
each connecting duct into two control volumes. The 
two control volume representation for the RE-1000 
eliminated the need to compute effective flow areas 
for those connecting ducts whose flow areas are 
nonuniform. The connecting duct model assumes 
adiabatic conditions. 


The regenerator heat-transfer and friction- 
factor correlations were revised. Earlier corre- 
lations were for wire screen regenerators. The 
new correlations are for Metex, knitted-wire 
regenerators. The heat-transfer correlation^ and 
the friction-factor correlation^ were based on 
experimental steady-state flow data. The correl- 
ations are 

0.588k p 5,66 Re 0, 76 

U ' 


The last modification made to the code was on 
the displacer gas spring model. An option was 
added to run the code with the displacer gas 
spring represented with either a control volume or 
a spring constant and damping factor. Previously 
only the control volume representation was used. 
Great difficulty was experienced in matching the 
predicted piston-displacer phase angle with the 
experimental data using the control volume gas 
spring representation. This difficulty seems to 
indicate that the displacer gas spring model is 
inadequate. This inadequacy might be due to the 
improper modeling of the gas spring hysteresis 
loss (no phase lag between heat transfer and aT 
is accounted for). Therefore, the spring constant/ 
damping factor option was chosen for the compar- 
isons made in this report. The spring constant 
and damping factor for the displacer gas spring 
are now calculated by the code during the con- 
strained mode of operation and are then used as 
inputs for the unconstrained mode. The spring 
constant is used in the calculation of the instan- 
taneous pressure of the gas spring: 


K . 

p h = p r + -r 1 x . 
d B A dr d 


The average pressure in the displacer gas spring 
is assumed to be equal to the average pressure in 
the bounce space. This assumption seemed reason- 
able since the displacer gas spring volume and 
the bounce space volume are in direct communica- 
tion every time the centering ports open. A damp- 
ing factor is used in the calculation of the net 
force acting on the displacer: 


F d ” ( P E P C^ A d + (P C V A dr + F B 


C, *, 
ds d 


Calibration Procedure 


The free-piston code can simulate three dif- 
ferent modes of operation: (1) constrained mode, 

(2) semiconstrained mode, and (3) unconstrained 
mode. In the constrained mode, the motions of 
both the piston and displacer are specified. This 
mode of operation is useful in calibrating the 
engine thermodynamics. In the semiconstrained 
mode, the motion of the piston is specified while 
the displacer motion is calculated. This mode of 
operation is useful in calibrating the displacer 
dynamics. In the unconstrained mode, the motions 
of both the piston and the displacer are calcu- 
lated. This mode of operation is useful in cali- 
brating the piston dynamics. 

Calibration Parameters 


The free-piston Stirling engine code had to 
be calibrated for the RE-1000 engine before com- 
parisons between predicted and experimental engine 
performance could be made. The code was calibrated 
by adjusting nine model parameters which are dif- 
ficult to specify accurately from geometrical 
data. These nine parameters are as follows: 


F = 12.27Re"°‘ 73 Re < 20 

F = 4.04Re -0 ’ 37 20 < Re <85 

F = 2.19Re -0 * 23 Re > 85 


FMULT pressure drop multiplication factor 

RHCFAC regenerator heat-transfer multiplication 
factor 

FHEXP expansion space cylinder heat-transfer 
multiplication factor 

FHCOMP compression space cylinder heat-transfer 
multiplication factor 
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CPTDO displacer centering port flow coefficient 

FDBLK displacer leakage multiplication factor 

CPTPO piston centering port flow coefficient 

FLEAK piston leakage multiplication factor 

CV2LD dashpot load coefficient 

Note that some of the parameters are labeled factor 
while other parameters are labeled coefficient. 

The calibration parameter list contains two dif- 
ferent subsets of parameters. The term factor 
refers to the subset of dimensionless calibration 
parameters while the term coefficient refers to 
the subset of calibration parameters that have 
physical properties. Calibration parameters are 
defined as the set of multiplication factors and 
coefficients required to adjust predicted pressure 
drops, heat transfer, and gas flow rates. 


displacer bore to the centering port in the dis- 
placer rod. For the piston, leakage occurs between 
the compression space and the bounce space. The 
major leak path in this case is the one between 
the centering port in the piston to the centering 
port in the cylinder wall. 

The dashpot load coefficient (CV2LD) is used 
to vary the piston stroke. Typically, CV2LD is 
varied until the code predicts the desired piston 
stroke. 

Initially, all of the calibration parameters 
are set at their nominal values. Running the code 
with the nominal calibration parameters represents 
a code without engine-specific factors. The 
nominal values are 


The pressure drop multiplication factor 
(FMULT) is applied to the pressure drops calcula- 
ted in the cooler, regenerator, and heater. The 
friction factor correlations assume steady flow 
conditions. Oscillating flow, which is the actual 
flow condition in the Stirling engine, is currently 
poorly understood; thus, effects of oscillating 
flow are not included in the friction-f actor (or 
heat-transfer) correlations. 


FMULT *1.0 
RHCFAC *1.0 
FHEXP = 1.0 
FHCOMP -1.0 
CPTDO = l.OxO -3 lb m /s 
FDBLK -1.0 
CPTPO = l.OxlO" 3 lb m /s 
FLEAK =1.0 
CV2LD = 2.25xl0~ 2 lbf-s 2 /in. (at engine design 
conditions) 


The regenerator heat-transfer factor (RHCFAC) 
is applied to the regenerator heat-transfer coef- 
ficient. The regenerator heat-transfer correla- 
tion also assumes steady flow conditions. RHCFAC 
was added to the code when the stacked wire screen 
heat-transfer correlation was being used to model 
the regenerator. Now that the code contains a 
correlation for Metex regenerators, RHCFAC is 
probably of less importance. 

The expansion and compression space cylinder 
heat-transfer factors (FHEXP and FHCOMP) are 
applied to their respective heat-transfer coeffic- 
ients. Modifications of these factors adjust the 
cylinder hysteresis (or cyclic) losses. There is 
a transient heat-transfer process between the work- 
ing space gas and the cylinder wall as the gas 
temperature is lowered and raised during the 
expansion and compression process. The net heat 
flow across this alternating temperature differ- 
ence gives rise to an irreversible loss known as 
the hysteresis loss. 

The centering port flow coefficients (CPTDO 
and CPTPO) are used to regulate the amount of 
centering port flow. They are used when the cen- 
tering ports are either completely or partially 
open. The purpose of the centering ports is to 
insure that the midpoints of the strokes of the 
displacer and piston stay at some fixed location 
in the engine. The centering ports are required 
since the leakage past most of the seals in the 
engine is not symmetrical. In the configuration 
of the RE-1000, the displacer will drift toward 
the compression space while the power piston will 
generally tend to drift into the working space. 

The leakage factors (FDBLK and FLEAK) are 
used to control the displacer and piston leakages. 
They are used whenever the centering ports are 
closed. Several leak paths exist for both the 
displacer and piston. The free-piston code cur- 
rently simulates the major leak paths. For the 
displacer, leakage occurs between the displacer 
gas spring and the bounce space. The major leak 
path is the one between the centering port of the 


The nominal value for the load coefficient CV2LD 
is calculated using the equation 


CV2LD = 


o P_JL sin e 
d cap 

32 ir 2 V 2 

U> A 


This equation gives the value of the velocity- 
squared load coefficient required to produce a 
sinusoidal piston motion with frequency a> and 
amplitude X p when the compression space pres- 
sure wave is sinusoidal with amplitude P ca 
and phase e relative to the piston motion. 


Constrained Mode Calibration 


The experimental case selected for the ini- 
tial step of the code calibration was run 1010 of 
the Lewis RE-1000 sensitivity test data. 3 This 
experimental data point was taken with the engine 
operating at its design conditions. The operating 
conditions were as follows: piston stroke^ 2.60 cm; 

heater-tube outside-wall temperature x 600 C; 
cooling-water inlet temperature, 25 C; mean 
working-space pressure, 7.0 MPa (1015 psi); work- 
ing fluid, helium. The calibration process began 
with the model operating in the constrained mode. 

The constrained mode was used to calibrate the 
code for the RE-1000' s thermodynamics. In this 
mode of operation the motions of the piston and 
displacer are specified as Fourier series. Only 
the fundamental terms (first harmonics) of the 
measured motions were used as inputs. The ampli- 
tudes of the higher order terms were found to be 
relatively small compared with the amplitudes of 
the fundamental terms which indicates that the 
piston and displacer motions were almost sinu- 
soidal. As a result, the difference between the 
code predictions using only the fundamental terms 
of the Fourier series and the code predictions 
using the first three terms of the Fourier series 
is insignificant (Figs. 5 and 6). This was found 
to be true for the RE-1000 with the dashpot load 
but may not be true for a different type of load 
(i.e., a hydraulic pump load). The inputs for the 
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piston and displacer in the form of Fourier series 
were as follows (amplitudes are in centimeters): 

x * -0.003 + 1.300 sin u>t 
P 

x d = 0.470 + 1.170 sin(ut + 57.5°) 

The factors FHCOMP, FHEXP, FMULT, and RHCFAC 
were available for adjustment in the constrained 
mode. These parameters were chosen since they all 
directly affect the engine's thermodynamics. FHCOMP 
and FHEXP were adjusted to obtain the correct com- 
pression space pressure wave amplitude and phase 
angle. FMULT was adjusted until good agreement 
was reached between predicted and measured pres- 
sure drop amplitudes across the cooler, regener- 
ator, and displacer. The experimental pressure 
drop amplitude measurements were fairly accurate 
as indicated by Schreiber et al . A great amount 
of work was done on the measurement techniques to 
insure accuracy. RHCFAC was left at its nominal 
value since the predicted indicated efficiency was 
already within the experimental error band. Once 
good agreement was reached between measured and 
predicted pressure drops, pressure amplitudes, and 
phase angles, the estimated cylinder wall tempera- 
tures were varied until good agreement was reached 
between the measured and predicted expansion and 
compression space gas temperatures. 

The displacer spring constant and damping 
factor required to achieve the correct displacer 
motion were then calculated by the code from a 
phasor diagram (Fig. 7) for use in the un con- r 
str ained mode pf operation. The vectors PcAdr 
and (Pc - PE) A d were calculated by Fourier 
analysis of the simulated pressure dr op and com- 
pression space pressure. The vector M^d was 
obtained from the differentiation of the displacer 
amplitude vector X^. The gas spring force 
required to produce the simulated displacer dynam- 
ics can be calculated from the following vector 
operation: 


F s = M d*d ~ t p C A dr + ( p C “ p F)A d ] 

The vector F s can then be resolved into the 
damping and spring components so that 


p s - C d s x d + ^ds x d 

The value of the displacer spring constant as 
calculated by the code was checked against the 
adiabatic spring constant as calculated using the 
equation 



This was done to insure that the spring constants 
calculated by the code were realistic in that they 
cannot be greater than their adiabatic limits. 

Semiconstrained Mode Calibration 

In this mode of operation, the motion of the 
piston is specified (constrained) while the dis- 
placer motion is calculated (unconstrained). This 
mode of operation aids in the calibration of the 
displacer dynamics when a control volume repre- 
sentation is used for the displacer gas spring 
model. But since a spring constant/damping factor 


representation was used, the semiconstrained mode 
of operation was not necessary for the calibration 
process documented in this report. 

Unconstrained Mode Calibration 

In the unconstrained mode of operation, the 
motions of the piston and displacer are calculated. 
CV2LD, FLEAK, and CPTPO are available for adjusting. 
These parameters are primarily used to adjust the 
piston dynamics. The parameters FDBLK and CPTDO, 
which are used to adjust the displacer dynamics, 
were left equal to their nominal values since the 
spring constant/damping factor representation was 
being used for the displacer gas spring. CV2LD 
was adjusted until the predicted and measured 
piston strokes matched. 

If the predicted and measured displacer 
strokes did not match, the factor FMULT (which 
was previously adjusted in the constrained cali- 
bration) had to be readjusted. Changing the value 
of FMULT changes the constrained calibration which 
means that the calibration process must be 
repeated. 

Several iterations between constrained and 
unconstrained modes were required to complete the 
calibration for the design point, run 1010. 

Table 1 compares the experimental data for the 
design point to the final predicted performance 
parameters for both the constrained and uncon- 
strained simulations. A flow chart summarizing 
the calibration process is shown in Fig. 8. 

After the calibration was completed for the 
design point, predictions were made for off-design 
conditions and were compared with experimental 
data. Each off-design point had only one of the 
operating conditions varied from the design con- 
ditions. The experimental runs used for the off- 
design conditions were as follows: (1) run 1006 

(piston stroke, 1.80 cm), (2) run 1012 (piston 
stroke, 3.00 cm), (3) run 1030 o (heater-tube 
outside-wall temperature, 450 °C), (4) run 1079 
(mean working-space pressure, 4.0 MPa), and (5J 
run 1200 (cooling-water inlet temperature, 55 C). 
These runs represented the extremes of the operat- 
ing conditions. The calibration parameters were 
then adjusted to improve the agreement over the 
entire range of operating conditions. The major 
calibration effort concerned the reference case 
or design point. Minor adjustments were then 
made for the five off-design points to improve 
the overall agreement. 

Results and Discussion 
Final Set of Calibration Factors 

The free-piston Stirling engine code was 
calibrated against a matrix of six experimental 
data points which represented the design point 
and off design points (extremes) of engine operat- 
ing conditions. The set of calibration param- 
eters which resulted in the best agreement between 
the predicted and measured data over the entire 
matrix of calibration runs was: 

FMULT = 1.2 
*RHCFAC =1.0 

FHEXP = 230.0 

FHCOMP = 50.0 

*CPTD0 = 1.0x10 J lb m /s ypsi 
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*FDBLK 1.0 , , 

*CPTP0 = l.OxlCT 3 lb m /s VPSi 
FLEAK = 1.2 

CV2LD = 2.24xl0” 2 lbf-s 2 /in. (Value of 
CV2LD for run 1010 only) 

where * indicates factors and coefficients still 
at their nominal values. 

Most of the final calibration factors and 
coefficients were very close to the nominal ones. 
Although the final values of FHEXP and FHCOMP appear 
quite large, when they are compared with the values 
required to isothermal ize the cylinder, the final 
values are relatively low. Figure 9 shows the 
expansion and compression space gas temperature 
ranges plotted as a function of the calibration 
factors FHEXP and FHCOMP. For nonisothermal 
cylinders, there is a large cyclic gas temperature 
variation in the expansion and compression spaces. 
The gas temperature range indicated in Fig. 9 is 
equal to twice the amplitude of this cyclic vari- 
ation. Therefore, a zero range indicates nearly 
isothermal conditions. Both FHEXP and FHCOMP were 
varied simultaneously over a wide range of values. 
The value of FHEXP required to isothermal ize the 
expansion space was found to be ~20 000 while the 
value of FHCOMP required to isothermal ize the com- 
pression space was found to be ~100 000. One pos- 
sible reason for the big difference between these 
two values is that in the RE-1000 the average sur- 
face area for heat transfer in the expansion space 
and the displacer appendix gap is almost four 
times larger than the average surface area for 
heat transfer in the compression space. It should 
be noted that for the calibration values of FHEXP 
and FHCOMP selected, the average heat-transfer 
coefficient predicted in the expansion and compres- 
sion space o at engine design conditions was 7.6 and 
6.8 W/cirr-°C, respectively. For the isothermal o 
case, these values would be 527 and 1300 W/cm 2 -°C, 
respectively. 

Additional Runs 

Four additional data points, runs 1017, 1024, 
1070, and 1121, were used as intermediate points 
for the final comparisons. 3 Runs 1017 and 1024 
were o for mean heater-tube temperatures of 550 and 
500 °C with the design operating conditions of 
piston stroke of 2.60 cm, cooling-water inlet tem- 
perature of 25 °C, and mean working-space pressure 
of 7.0 MPa. Experimental run 1070 was for an 
intermediate mean pressure of 5.5 MPa while run 
1121 was for an intermediate cooling-water inlet 
temperature of 40 °C with the other operating 
conditions at their design values. 

Comparison of Engine Dynamics and Thermodynamics 

The final set of calibration parameters was 
used to make comparisons of predicted piston and 
displacer strokes, phase angles, and frequency; 
gas pressure amplitude; and qas temperatures in 
the expansion and compression spaces against 
measured data. Comparisons of predicted and 
measured data are shown as a function of piston 
stroke in Figs. 10 through 16. The piston stroke 
of the engine was adjusted by varying the engine 
load; the piston stroke in the simulation was 
adjusted by varying the dashpot load coefficient 
(CV2LD). The comparisons are shown for a range 
of piston strokes with the heater-tube temperature, 
cooling-water inlet temperature, and mean working 


space pressure set at their design values. Error 
bands have been placed on the experimental data. 

The error bands indicate the measurement error 
associated with each reading of data. 3 The 
ultimate goal of the calibration process was for 
the code predictions to fall within the error 
bands. 

Figure 10 shows the predictions of the piston 
and displacer strokes to be in good agreement over 
most of the tested load range. Figure 11 shows 
the predicted piston-displacer ghase angle within 
the experimental error band (^1°) over the entire 
load range. This good agreement can be attributed 
to the use of the spring constant/damping factor 
representation of the displacer gas spring. It 
should be noted that the actual piston-displacer 
phase angle at the engine design conditions was 
higher than the design phase angle of 45°. This 
difference was probably due to a slightly higher 
displacer damping force than assumed during the 
design process. A phasor diagram of the forces 
acting on the displacer was helpful in reaching 
this conclusion (Fig. 7). The magnitude of the 
displacer damping force is difficult to predict. 

The phasor diagram shows that the displacer damp- 
ing force has a direct affect on the piston- 
displacer phase angle. Figure 12 shows the 
predicted piston-pressure phase angle falls within 
the experimental data's error band (*1°) over most 
of the tested load range. The good agreement 
between the predicted and measured phase angle is 
a result of increased cylinder heat transfer by 
means of the calibration factors FHEXP and FHCOMP. 
These two factors also had a positive effect on 
the compression-space pressure amplitude. The 
compression-space pressure amplitude has always 
been overpredicted in the past. Increasing the 
cylinder heat transfer brought the predicted 
amplitude within the error band (±20 kPa) of the 
experimental data. Figure 13 shows that the pre- 
dicted engine frequency increases slightly with 
stroke while the measured frequency is almost 
constant. Figures 14 and 15 show the expansion 
and compression space gas temperatures compared 
with the experimental data. The good agreement 
shown for the expansion space can be attributed 
to the better estimates of the cylinder-wall tem- 
peratures. The good agreement shown for the com- 
pression space is a result of the new cooler model 
and the better estimates of the cylinder-wall tem- 
peratures. Figure 16 shows the predicted 
compression-space pressure amplitude to fall 
within the experimental data's error band (±20 kPa) 
over most of the tested load range. 

Comparison of Power and Efficiency 

Figures 17 to 24 show comparisons for power 
and efficiency between the engine test data and 
the code predictions. These comparisons were made 
for varying piston stroke, mean working-space 
pressure, mean heater-tube temperature, and cooling- 
water inlet temperature. Fiqure 17 shows the pre- 
dicted power to fall within the experimental error 
band (±40 W) over most of the range of piston 
strokes. Figure 18 shows the predicted efficiency 
to fall outside the error band (±1.3 percent) at 
the lower piston strokes. The experimental data 
show a significant drop in the engine efficiency 
at the high piston stroke. One possible reason 
for this might be due to a higher heat-transfer 
loss out of the expansion space at the larger 
strokes. As the piston stroke is increased, the 
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displacer stroke is also increased. This allows 
the expansion-space gas to come into contact with 
a cooler portion of the cylinder wall. The model 
does not simulate the drop in efficiency because 
the expansion-space wall temperature is assumed to 
be constant over its entire length. Comparisons 
of measured and predicted power and efficiency are 
shown as a function of mean working-space pressure 
in Figs. 19 and 20. The comparisons are shown for 
a range of mean pressures with the piston stroke, 
heater-tube outside-wall temperature, and cooling- 
water inlet temperature set at their design values. 
Figure 19 shows the predicted power to be slightly 
above the experimental data over most of the pres- 
sure range tested. Figure 20 shows the predicted 
efficiency to fall within the experimental error 
band (±1.3 percent) over the entire range of pres- 
sures tested although the shapes of the curves 
are significantly different. 

Comparisons of measured and predicted power 
and efficiency are shown as a function of mean 
heater-tube outside-wall temperature in Figs. 21 
and 22. The mean heater-tube outside-wall temper- 
ature was experimentally measured by taking the 
average of twelve heater-tube temperatures; six of 
which were measured at the quarter-length point at 
the expansion space end of the heater tube, while 
the other six were measured at the quarter-length 
point of the heater tube near the regenerator end 
of the tube. For these comparisons the piston 
stroke, mean working-space pressure, and cooling- 
water inlet temperature remain at their design 
values. Figure 21 shows the predicted power to 
fall within the experimental error band (±40 W) 
over the entire range of heater-tube temperatures. 
Figure 22 shows the predicted efficiency to fall 
below the experimental data over most of the range 
of heater tube temperatures although the shapes 
of the curves are similar. 

Comparisons of measured and predicted power 
and efficiency are shown as a function of cooling- 
water inlet temperature in Figs. 23 and 24. For 
these comparisons the piston stroke, mean working- 
space pressure, and mean heater-tube outside-wall 
temperature remain at their design values. 

Figure 23 shows the predicted power to fall 
slightly above the experimental data over most of 
the range of cooling-water inlet temperatures. 
Figure 24 shows the predicted efficiency to fall 
within the experimental error band (±1.3 percent) 
over the entire range of cooling-water inlet tem- 
peratures. 

Concluding Remarks 

Many modifications have been made to the NASA 
Lewis free-piston Stirling engine code. The 
recent improvements have had a positive affect on 
the code predictions. Further development of the 
free-piston code should bring the predicted and 
measured data into even better agreement. The 
ultimate goal is to have a code that does not 
require any engine-specific calibration parameters. 
This would allow the code to be used to model any 
free-piston Stirling engine without having to go 
through the calibration process. More importantly 
the code could be used in conjunction with design 
codes to confidently design future free-piston 
Stirling engines. The results of the work done 
and documented in this report indicate several 
areas that require improvement: 


1. New heat-transfer correlations are required 
for the expansion and compression spaces. This 
heat transfer is not well understood and the 
effect appears to be significant. 

2. The friction-factor and heat-transfer 
correlations for the three Stirling cycle heat 
exchangers assume steady-state conditions. The 
effects of oscillating flow and pressure level 
need to be better understood. 

3. The leakage and centering port flow cal- 
culations should be verified. Currently, no 
experimental data of this type exist for the 
verification. 

4. More information is required on gas 
springs in terms of irreversible heat-transfer 
(hysteresis) loss. 

Current plans are to continue the validation 
of the free-piston model for the RE-1000 engine 
with a hydraulic output load. A subroutine for 
simulating the hydraulic load has been added to 
the free-piston code. The hydraulic output sub- 
routine needs further development. The RE-1000 
free-piston Stirling engine hydraulic output 
system is currently being debugged and readied 
for testing at NASA Lewis. 

Other code validation work is also proceeding 
at NASA Lewis with the Stirling engine kinematic 
code. The predictions of this code with no 
engine-specific calibration parameters are being 
compared to test data from a variety of Stirling 
engines. RE-1000 test data and predictions of the 
free-piston code as described in this report but 
without the engine-specific calibration parameters 
will be input to this validation effort. 
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TABLE 1. Comparison of RE-1000 Predicated and Experimental Performance at Design Operating Conditions (Run 1010) 


Run 

description 

Frequency, 

Hz 

Heat in, 
kW 

Heat out, 
kW 

Piston power 

Dash- 

pot, 

kW 

Indicated 

efficiency, 

percent 

Cycle 

energy 

balance 

error 

Piston 

energy 

balance 

error 

Indicated PV, 
kW 

Phasor, 

kW 

Experiment 

Predicted: 

30.1 

3.643 

2.736 

0.939 

1.021 

0.866 

25.8 

0.88 

-7.77 

Constrained 

30.1 

3.921 

3.036 

.939 

.944 

.974 

23.9 

1.38 

3.73 

Unconstrained 

30.3 

3.939 

2.887 

.969 

.951 

.994 

24.6 

-2.12 

2.58 


P A V cos \|> 

Experimental indicated power = ca 


Predicted indicated power is determined by integrating P-V diagram 
Phasor diagram power = na,p ca ApXp sin e 

* * • . . . Indicated power 

Indicated efficiency = Heat 

~ , -_ A Indicated power + heat out . 

Cycle energy balance error = 100 1 

Power piston energy balance error = 100 [ DaShp ° * d i ' c ate^ power** ' ° ~ J ] 
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OFF-DESIGN DESIGN OFF-DESIGN 

PISTON STROKE, 1.8 CM PRESSURE, 4.0 MPa 

ENGINE OPERATING CONDITIONS 

FIG. 5 COMPARISON OF INDICATED POWER PREDICTIONS USING 
ONE AND THREE FOURIER SERIES TERMS FOR PISTON AND DIS- 
PLACER. 



FIG. 6 COMPARISON OF INDICATED EFFICIENCY PREDICTIONS 
USING ONE AND THREE FOURIER SERIES TERMS FOR PISTON 
AND DISPLACER. 
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FIG. 7 DISPLACER PHASER DIAGRAM 
OF RUN 1010 (REFERENCE CASE). 
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- LEWIS CONSTRAINED MOTION SIMULATION 
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FIG. 8 FLOW CHART OF CALIBRATION PROCESS. 
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DISPLACER STROKE. CM GAS TEMPERATURE RANGE 



CYLINDER HEAT-TRANSFER FACTORS FHEXP AND FHCOMP 

FIG. 9 EXPANSION AND COMPRESSION SPACE GAS TEMPERA- 
TURE RANGES VERSUS FHEXP AND FHCOMP (RANGE = 2 x 
AMPLITUDE). 



FIG. 10 DISPLACER STROKE VERSUS PISTON STROKE. 
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FIG. 16 COMPRESSION-SPACE PRESSURE AMPLITUDE VERSUS 
PISTON STROKE. 
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INDICATED EFFICIENCY, PERCENT j INDICATED POWER, W INDICATED EFFICIENCY, percent 



MEAN HEATER-TUBE OUTSIDE-WALL TEMPERATURE, °C 

FIG. 22 INDICATED EFFICIENCY VERSUS MEAN HEATER-TUBE 
TEMPERATURE. 



FIG. 23 INDICATED POWER VERSUS COOLING-WATER INLET 
TEMPERATURE. 



COOLING-WATER INLET TEMPEATURE, °C 

FIG. 24 INDICATED EFFICIENCY VERSUS COOLING-WATER INLET 
TEMPERATURE. 
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